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A simple, causal, and analytic model of final state rescattering is used to describe all available data 
on the exotic resonances Zc(3900) and Zc(4025). The model provides a compelling and accurate 
representation of experiment with no need for poles in the scattering matrix. 
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I. INTRODUCTION 

The recent confirmation of the Z(4475) resonance by LHCb[T] and the discovery of charged charmonium states 
Zc(3900), Zc( 4025) and bottomonium states, Z;,(10610), Zf,(10650)[2], point to a possible extensive exotic hadronic 
spectrum. These discoveries have led to a host of interpretations of these states as loosely bound molecules [3], 
tetraquarks[4], hadrocharmonium[5], or hybrids [6]. This paper explores the simplest possible explanation of the 
Zc(3900) and Zc(4025), namely that the enhancements associated with these states are associated with threshold or 
cusp effects. 

The Zc(3900) was discovered by the BESIII and Belle collaborationsin e+e” —>■ F(4260) —>■ J/V'tt+tt” in the 
charged mode Zc —>■ The reported mass and width are M = 3899.0 ± 3.6 ± 4.9 MeV and T = 46 ± 10 ± 20 

MeV. There is strong evidence that the quantum numbers of the state are = I"*". The Z’c(3900) was also observed 
as a threshold enhancement in the reaction e“'"e“ —>■ ttDD*^, where its mass and width were determined to be 
3883.9 ± 1.5 ± 4.2 MeV and 24.8 ± 3.3 ± 11.0 MeV respectively. 

The Zc(4025) was observed by BESIII in e+e” —)■ D*D*tt at ^/s = 4.26 GeV[9] and in e+e” —>■ /icTttt at a variety of 
energies from = 3.90 to 4.42 GeVlTU]. Its mass was determined to be 4026.3 ±2.6 ±3.7 MeV and 4022.9 ±0.8 ±2.7 
MeV in the respective experiments, while the measured widths were 24.8 ± 5.6 ± 7.7 MeV and 7.9 ± 2.7 ± 2.6 MeV 
respectively. 

Gaining an understanding of this spectrum is clearly relevant to advancing qualitative (and possibly quantitative) 
understanding of quantum field theories in their nonperturbative guises. In this regard, it is important to develop 
a sufficiently robust phenomenology to be able to distinguish perturbative and nonperturbative explanations of ex¬ 
perimental data. Eor example, it has been noted that cusp effects due to coupled channel thresholds can lead to 
enhancements in rates that mimic S-matrix poles[TTJ [12]. This point was recently expanded in Ref. |13| . where it 
was shown that a simple model of a four-point function that incorporates a virtual coupled channel [DD* and D*D* 
or BB* and B*B*) can provide a unified and quantitatively accurate description of the four exotic states Zc(3900), 
Ze(4025), Zb(10610), and ^6(10650). 

The conclusions of Refs. mma were subsequently challenged by Guo et aZ.[T3], who argued that explaining the 
processes V(4260) —>■ ttDD* and V(4260) —)■ imJ/ip required nonperturbative interactions, which naturally generate 
S-matrix poles and obviate the cusp explanation of the charged charmonium states. This conclusion was based on 
a model that assumed pointlike coupling between relevant hadrons (tt, V(4260), D, D*, and J/tp). The threshold 
enhancement in V —>■ ttDD* was then generated via a loop with an intermediate DD* state that was regulated 
with a phenomenological exponential cutoff function. The resulting coupling constants were too small to provide 
agreement with Y —>■ TnvJ/4> data unless the loop diagrams were iterated, which led the authors to conclude that the 
Zc enhancements are true resonances. 

The purpose of this paper is to examine the robustness of the conclusions of Guo et al. by constructing an explicit 
model for the interactions of the relevant hadrons. This will also permit checking the conclusions of the model 
presented in Ref. m- The proposed model is both simple and plausible. Furthermore, its predictions follow from 
the Schrodinger equation, thereby avoiding the criticism of Ref. |15| . Model predictions are in good agreement with 
the available charmonium data and provide strong evidence that dynamically generated poles are not necessary to 
explain the novel charged charmonium states. 


II. COUPLED CHANNEL MODEL 

The proposed model will employ nonrelativistic dynamics and separable hadronic vertices. It is a simple matter to 
relax these choices; however, doing so will not affect results in a qualitative way, and we Hnd it preferable to compute 
with the simplest model possible to illustrate the relevant effects. The model vertices are of the form 
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AB-.CD-^ gAB-.CD ■ k‘^AB Pab) ' k^co ^^V^-klo!P cd)^ ( 1 ) 

where g is the coupling for the relevant channel, /3 is a hadronic scale of order Aqcd > and £ is an integer specifying an 
angular momentum factor. With the exception of £h^-K = Ij all other channels are dominated by S-wave scattering, 
hence £ = 0. The momenta, k, are given by 


kAB = HsAB]m\, m%)/{4:SAB) ( 2 ) 

where A is the Kallen function and sab = {pa+Pb)^ ■ This model will be used for the channels Ytt : DD*, Ytt : D*D*, 
J/ijj-K : DD*^ J/ijjTr : D*D*, heir : DD*, and heir : D*D*. Appropriate charge conjugation is included where required. 
We remark that these vertices are often represented as t-channel D or D* exchange diagrams. 

Rescattering will be accommodated with similar vertices for the channels DD* : DD* and D*D* : D*D*. Cross 
channel scattering DD* -A D*D* is not considered in the following. In other models, these vertices are commonly 
described by pion-exchange diagrams. It will be assumed that a direct coupling, Ytt : J/ipTr, is negligible; hence this 
process will be generated by iterations of diagrams of the type shown in Fig. [fright). 

Exponential form factors are adopted in the model vertices because hadronic constituents are not resolved at the low 
energies being considered here. We regard this assumption as more appropriate for the study of production than 
constant couplings. In general, each of the exponentials can assume a different scale, Pab- However, for simplicity 
we shall assume that the fis take on a single value for almost all of the form factors (the exception will be that for 
the DD* form factor, as described below). 


W\AA^ 






FIG. 1: Leading order diagrams for Y -A ttD*D{D*) (left) and Y -A TnTj/'tp{hc) (right). 


Fig.[T] shows the leading order diagrams for the processes considered here. All of the results presented below will 
be unitarized by summing the bubble diagrams that contribute in s-channel. The effect of this unitarization can be 
strong: sufficiently attractive couplings lead to a DD* or D*D* bound state and give a sharp peak below threshold 
in the appropriate channel. Similarly, a repulsive coupling smears out the threshold cusp and pushes it to higher 
invariant mass. It is thus possible to determine the size of these effects by comparison to experiment. 


III. CHARMONIUM CUSP STATES 


The model scales Pytt, Pdda Pd*d* and the couplings gy-K-.DD*, PYj-.d^da goD^-DD* gD*D*-.D*D* will be fixed 
by comparison to the threshold behavior of the processes Y -A ttDD* and Y -A ttD*D*. Since the purpose of this 
computation is not a detailed analysis of experimental results, a fit to the Dalitz plot densities will not be attempted, 
except where required. 


A. Y ^ ttD*D* 

Fig. I displays the projection of the Y -A ttD*D* Dalitz plot onto the m{D*D*) axis. The overall scale is set by 
the coupling pytt-.d-^d*, which is arbitrary since cross sections were not measured. The solid line shows the result of 
a fit that evidently describes the data quite well. The scales chosen were (iyir = /3d*d* = 0.3 GeV. This value is 
strongly selected by the data since deviations lead to quite different shapes. Similarly, D*D* rescattering also affects 
the predicted rate and the result 


gD-D‘:D‘D* ~ 0 


( 3 ) 
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is strongly preferred. Thus the reaction V —> ttZ)*!)* implies that no dynamical D*D* resonances will be generated 
in any of the following work. We also observe that there is no need for a constant (coupling) term in the amplitude, 
in opposition to what was assumed in Ref. [14] for Y —> ttDD* . 



FIG. 2: e''"e {s/s = 4.26 GeV) —>■ ttD*D*. Solid line: model fit. Data from Ref. |9]. 


B. y tvDD* 

Data for Y —>■ tyDD* are presented in Fig. In attempting to fit this, we note that only a pion and a single D 
meson were reconstructed during data-taking - other particles were inferred [7] . Thus it is likely that an incoherent 
background exists in this data and we therefore choose to model this reaction by incorporating such a background as 
a constant. The scale Pytt has already been fixed to 0.3 GeV; this leaves (Sdd* and the coupling qy-k-.dd* (which sets 
the normalization) and gDD*-.DD* to determine from the data. 

A rough fit (no systematic fitting was made in this work as it was thought to be antithetical to the approach) 
yielded I5dd* = 0.2 GeV (shown as a solid line in Fig. |^. The fit is gratifyingly faithful to the data. 

Once again, the data present no evidence for strong rescattering, and we determine goD^-.DD* ~ 0, implying that 
no DD* resonances can be dynamically generated. These conclusions disagree with those of Guo et al., who assume 
that a constant gYir-.DD* vertex gives rise to the events seen at large DD* invariant mass. As we have seen, this 
is at odds with the data for Y ■nD*D* and is not necessary to describe Y —>■ ttDD*. Furthermore, preliminary 
results from BESIII do indeed indicate that the background is strongly suppressed when all three final particles are 
reconstructed [IB] . providing further support for the model presented here. 


C. Y^TT+lY-J/lp 

With all the relevant scales and the D^*^D* couplings fixed, it is possible to make a prediction for Y TnrJ/ij;. 
Unfortunately, this process is complicated by strong final state interactions in the tttt subsystem, shown in the left panel 
of Fig. 1^ Thus a reasonably faithful reproduction of the data requires incorporating tttt dynamics in the amplitude 
model. Because pion dynamics is not the thrust of the current investigation, a pair of Breit-Wigner amplitudes at 
0.35 GeV and 0.88 GeV with widths of 68 and 290 MeV respectively were used to obtain an approximate fit to the 
TTTT spectrum (shown as a solid line in Fig. [^. A Flatte parameterization of the /o(980) was also attempted but this 
did not fit the data well. 

A constant coherent background is assumed and the assumption gj/^Tr-.oD* = ‘^gj/ip-R-.D* d* is made. With the model 
of TTTT dynamics in place, the amplitude is fixed (including the overall normalization, which is determined by the tttt 
distribution). The resulting prediction is shown as a solid line in the right panel of Fig. The overall quality of 
the prediction is quite good, although some strength is missing through the .Zc(3900) reflection near 3450 MeV. Note 
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FIG. 3; e^e (-^i = 4.26 GeV) —>■ -kDD*. Solid line: Pdd* = 0.2 GeV. Data from Ref. [S]. 




FIG. 4: e''"e {-s/s = 4.26 GeV) —> 'KiiJ/ip. Left panel: invariant tttt mass distribution. Right panel: invariant 'Kj/'ij) mass 
distribution. Filled squares: 'r~ J/ip\ open squares: Jjrp. Data from Ref. [7]. 


that there is a slight enhancement in the prediction that corresponds to the Zc(4025) cusp in D*D*. The relative 
size of this enhancement is controlled by the ratio 5 j/i/, 7 r:_D_D*/ 5 j/^ 7 r:D*D* ■ This is expected to be of order unity. As 
mentioned above, the ratio has been set to 2, but in principle could be larger, so that the D*D* cusp would be barely 
visible in the plot. 


D. e^e —>■ 7r"'"7r he 

The ^c( 4025) was observed in e+e” —>■ inrhc, which was measured by BESIII at 13 values of ^/s. The total event 
rate summed over all energies is reproduced as data points in Fig. The data were modelled by generating 13 Dalitz 
plots corresponding to the experimental values of y/s and summing these with a weight given by the reconstructed 
number of mesons for each energy. All couplings and scales are fixed as above, with the exception of which 
was set to 0.3 GeV for simplicity. 

In principle, tttt rescattering should affect the predicted distribution. However, we have found that employing the 
fit dynamics from Y —>■ tttt J/ip does not change the results much; thus a simple constant coherent background was 
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m{ii; hj (GeV) 


FIG. 5: e^e —>■ h^. Solid line: model fit. Data from Ref. |10 |. 


employed in the fit shown here. Finally, the relative couplings were fit to the data to obtain gh^-K-.D*D* ~ ‘^gh^-w-.DD* ■ 
Note that the analogous ratio is approximately inverted in the case of the J/'i/'TTTr final state. 

Once again, the agreement with data is quite satisfactory, especially considering the simplicity of the model. The 
peak near 4.0 GeV is a cusp effect due to D*D* rescattering; the peak near 3.9 GeV is due to DD* rescattering, and 
the peak near 3.75 GeV is a reflection of the Zc(4025) cusp. None of these features are associated with poles in the 
S-matrix. 


E. gN ^ fi J/‘ip Tv^ N 

Lin et al. have suggested that the coupling of the Zc(3900) to TTj/ij) can be exploited to search for it in 
photoproduction |17j . The idea is that the virtual photon converts to a J/ip via the vector meson dominance mech¬ 
anism, which interacts with a nucleon by pion exchange, creates an s-channel which finally decays to TrJ/'i/;. 
The cross section for yV ZN was estimated using a hadronic Lagrangian with dipole form factors. The Z^itJ j'lp 
coupling was taken from the measured width of the Z,. (assuming that it is saturated by this mode). The resulting 
cross section was predicted to peak at ^/s « 7 GeV with a readily observable rate. 

In spite of these expectations, a measurement of gN —>• gJ/i/jirN by the GOMPASS collaboration [18] found no 
evidence for the Zc(3900). This lack of evidence has a simple explanation in the present model. Dipole form factors 
should be replaced with the form factors of Eq. which are heavily suppressed by the large center of mass energy 
of the process. Of course, this observation has no bearing on whether the Zc is a dynamically generated resonance 
since the coupling to a resonance state could be similarly suppressed. Furthermore, the model vertices should not be 
trusted at very large center of mass energies since this is where the pomeron trajectory is expected to take over the 
dynamics. 


IV. CONCLUSIONS 

A detailed comparison to available Zc experimental data has been presented. The model successfully captures 
the features of all the data and indicates that there is no evidence for strong DD* or D*D* rescattering in this 
system. We note that this conclusion is supported by recent lattice gauge computations that report only weakly 
repulsive {D*D*)^ interactions in the = I"*" channel [TU]. Thus isovector rescattering is not sufficiently attractive 
to generate dynamical bound states (this is also supported by lattice computations [20]) and exotic resonances are 
not required to explain the data. However, cusp effects remain and can provide a qualitative, and even quantitative, 
explanation of all of the data. 
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This simple model is fully quantum mechanical and therefore removes concerns about causality [IS]. Furthermore, 
it appears to provide a better, and more complete, description of the available data than that of Ref. [Hi- 

Threshold enhancements and openings are generic features in hadronic systems and one must therefore be cautious 
in claiming bound states where such effects are known to operate. Certainly, any near-threshold enhancement can 
simply arise because hadrons are soft; thus the DD* and D*D* data are easily explained. Similarly, coupled-channel 
cusps should be regarded as the possible explanation for bumps seen in rescattering channels slightly above coupled 
channel thresholds. Furthermore, if the “widths” of these enhancements vary strongly (as they do for the Zcs) 
between threshold and rescattering processes, then this is an additional sign that non-resonant explanations must 
be considered. In particular, threshold bumps arise due to competing effects between form factors and phase space, 
whereas a rescattering enhancement width is mediated by form factors and a rescattering loop. Importantly, this 
implies that (cusp dominated) threshold bumps do not exhibit phase motion, while rescattering enhancements may 
have phase motion due to the associated bubble diagrams. 
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